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Abstract 

We show that the well-known divergence of the perturbative expansion of resummed results for 
processes such as deep-inelastic scattering and Drell-Yan in the soft limit can be treated by Borel 
resummation. The divergence in the Borel inversion can be removed by the inclusion of suitable 
higher twist terms. This provides us with an alternative to the standard 'minimal prescription' 
for the asymptotic summation of the perturbative expansion, and it gives us some handle on the 
role of higher twist corrections in the soft resummation region. 
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It has been known for some time [1] that the resummation of logarithmically enhanced 
contributions to the coefficients of the QCD perturbative expansion due to soft gluon radiation 
has the effect of rescaling the argument of the strong coupling constant: the hard perturbative 
scale is replaced by a relatively soft scale related to the radiation process. For example, for a 
physical process characterized by the hard scale Q 2 and a scaling variable < x < 1 near the 
x — 1 boundary of phase space (e.g. close to the production threshold for a given final state), 
the resummation of large logs of 1 — x effectively replaces the perturbative coupling a s (Q 2 ) with 
a s {Q 2 {l — x)). Similarly, in the resummation of soft px spectra [2] the argument of the strong 
coupling becomes p T , and so on. Therefore, the perturbative approach eventually fails in this 
'soft' kinematical region. This failure is understandable on physical grounds, because at the 
phase space boundary, i.e. as x — > 1, the center-of-mass energy is just sufficient to produce the 
given final state, so, for instance in the case of deep-inelastic scattering in this limit the process 
becomes elastic. 

In practice, this problem must be treated in some way in order to obtain phenomenological 
predictions: eventually, at some low scale A (the position of the Landau pole) the strong coupling 
blows up, so when 

A 2 

x = Xl = 1 ~q2 (!) 

resummed results become meaningless. The scale A is usually identified with Aq CD . However, 
rr-space resummed results are known to run into difficulties anyway, regardless of the size of 
the coupling constant, essentially because the resummation of leading (or next fc -to-leading) 
logarithmic contributions in x space does not respect momentum conservation [3] : this produces 
a spurious factorial divergence of resummed results when expanded at fixed perturbative order. 
Because of this factorial divergence, attempts to remove the problem of the Landau pole by 
cutting off the resummed x-space result at large x & xl [4] display a sizable dependence on 
the choice of cutoff. It has therefore been suggested [3] that it is in fact more advisable to 
consider resummed results in terms of the variable A" which is Mellin conjugate to x, namely, 
to consider the resummation of large In A" as A" — > oo. Indeed, it can be shown that to any 
logarithmic order [5] upon inverse Mellin transformation the resummation of next fc -to-leading 
In A" contributions provides resummation of next fe -to-leading ln(l — x) terms, up to subleading 
contributions. 

If the A^-space resummed result is expanded perturbatively in powers of a s (Q 2 ) and then 
Mellin transformed back to x order by order, one winds up with a divergent series: in other 
words, the A^-space resummed result cannot be obtained as the Mellin transform of a pertur- 
bative x-space calculation. However, it is possible to give a "minimal prescription" [3] for the 
reconstruction from the A^-space resummed result of an x-space result to which this divergent 
sum is asymptotic (in the sense of asymptotic series). The minimal prescription (MP) thus leads 
to a result which is well-defined and smooth for all < x < 1. It has been widely used for 
phenomenological applications. 

Here, we wish to reconsider this issue, partly motivated by the results of ref. [5], which give 
us full control on the relation between large x and large A" resummations, and by the observation 
that the fact that the minimal prescription is well defined for all x is in fact a mixed blessing, 
if one cannot control what is happening when x ~ xl, where the perturbative approach breaks 
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down. We will see that the divergence of the x-space perturbative result can be traced to 
subleading terms when the Mellin inversion of the iV-space result is performed to all logarithmic 
orders (but excluding subleading powers), and that it can be treated by Borel resummation, 
at the expense of including higher-twist contributions. The result which is obtained shares 
the pleasing features of the minimal prescription: it provides an asymptotic summation of the 
divergent perturbative expansion, and it is well-defined for all x. However, it differs from it, 
though this difference only becomes significant for x ~ xl i.e. around the Landau pole (where 
it is in fact closer to the truncated perturbative result). Also, it has a rather different physical 
interpretation: if suitable higher twist contributions are included, the perturbative expansion 
becomes convergent, and our result provides its sum. 

We start by recalling how the Landau pole appears and is treated with the minimal prescrip- 
tion. Consider a generic observable <j(Q 2 , x) and its Mellin transform 

a{Q\N)= f 1 dxx N - 1 a(Q 2 ,x), (2) 
Jo 

which maps the region of large x onto the region of large N. The resummation of logarithms of 
N can be performed [5] in terms of the physical anomalous dimension 

d In a{Q\N) 

lMQ),N)= — 2 . (3) 

For example, in the case of deep-inelastic scattering structure functions the resummed expression 
of 7(a s (Q 2 ), N) has the form [5-7] 



lMQ 2 ),N) = / ™ J29ka k s (Q 2 /n) + 0(N°), (4) 
Ji n f-r. 

k=l 



where gt are constants to be determined by matching with fixed-order calculations, and the 
neglected terms are either iV-independent, or suppressed by inverse powers of N for large N. 
Truncating the sum at k = 1,2,... corresponds to computing 7(a s (Q 2 ), N) to leading, next- 
to-leading, . . . logarithmic accuracy. 1 Starting from the resummed expression of , y(a s (Q 2 ), N), 
one can obtain the Mellin transform of the cross section, resummed to the same logarithmic 
accuracy. The physical quantity cr(Q 2 , x) may be obtained by inversion of the Mellin transform. 

Equation (jlj) shows explicitly that the resummed result depends on a s (Q 2 /N): the resum- 
mation has replaced the hard scale Q 2 with the softer scale Q 2 /N — in fact, in the soft limit 
the resummed result only depends on N through this rescaled coupling. As a consequence, 
r j{a s {Q 2 ), N) has a branch cut along the real positive axis for N > N^, where Nl is the location 
of the iV-space Landau pole 

O 2 

^ = £ ■ (5) 

This in particular implies that the inverse Mellin transform of the resummed anomalous di- 
mension eq. does not exist. Indeed, the Mellin transform of a function f(x) such that 



1 Beyond next-to-leading log, eq. only holds provided the cross section has a particular factorization prop- 
erty. This is immaterial for the ensuing discussion, where we will concentrate on the leading logarithmic case. 
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|/(x)| < Kx~ N ° for all x, with K and N real constants, is an analytic function of the complex 
variable N in the half-plane ReiV > N . Therefore, ^(a s (Q 2 ),N) eq. (j3J) cannot be the Mellin 
transform of any function. 

To see the problem more clearly, let us consider for definiteness the resummed expression of 
7(a; s ((5 2 ), N) to leading logarithmic accuracy, 



7llK(Q 2 ), N) = 9ij"^ u s (Q 2 /n) = In (l + (3 a s {Q 2 ) In i 

where we have consistently used the leading- log expression of a s : 

a s (Q 2 ) _ 33 -2n f 



12tt 



(6) 



(7) 



Because of the Landau singularity, 7ll(«s(Q 2 ), N) has a branch cut on the real positive axis for 



N > N L = e?o<>s(Q 2 ). (8) 

One may formally consider the term-by-term inverse Mellin transform of the expansion of 
1ll(®s(Q 2 ), N) in powers of ot s {Q 2 ). This gives 



PllMQ 2 ),x) 



lim 



A 

E 



•N+ioc 1 

$<x k s(Q 2 )-^ I- dNx- N \n k -- N>0. (9) 

2m JN-ioo N 



Each term of the series is now a well-defined inverse Mellin transform, but the series does not 
converge, so we cannot take the limit K —>■ oo. Indeed, if the series were convergent, one could 
interchange the sum over k and the integral over N in eq. 0, but the sum 



E 

k=i 



A\ k+1 1 



is only convergent for 



Poa s (Q 2 ) In - 



< 1, 



(10) 



while the integral in iV on the path Re iV = N involves values of N outside this range. 

However, as well known, if the Mellin transform in eq. Q is computed at the relevant (leading, 
next-to-leading...) logarithmic level the perturbative series converges. Indeed, considering again 
for definiteness the leading log case one has 



dNx- N \n k — 
2ni JN-ioo N 



In^Hl-x) 



x 



+ NLL, 



(12) 



where + denotes the standard prescription of Altarelli-Parisi evolution [8]. The series in eq. 
now is convergent for all x < x^. its sum is 



PhL x (, a s{Q 2 ),x) = -9i 



a s (Q 2 



1- xl + (3 a s (Q 2 )\n(l- x) 



-9i 



a s {Q 2 {l-x)) 
1 — x 



(13) 



which is singular at the Landau pole eq. This singularity sets the radius of convergence of 
the series. Similar arguments can be used to show [5] that if we start with the next fc -to-leading 
In i result and we perform the Mellin inversion to the same order we wind up with a series 
of terms of the form of eq. (fl~3"j) . but with higher powers of the coupling, up to a k {Q 2 {l — x)). 
Equation (J13|) shows explicitly that the resummation replaces the scale Q 2 with Q 2 (l — x). 
Notice that eq. (|13|) and its higher-order cognates provide us with a next fc -to-leading ln(l — x) 
resummation at the level of the physical anomalous dimension (or rather splitting function), 
rather than of the cross section, and it is free of the spurious factorial growth mentioned above. 
In fact, it can be shown that the factorial growth is a byproduct of the next fc -to- leading log 
truncation of the exponentiated result, and it disappears provided only the Mellin transform 
of the exponentiated result is determined including subleading logarithmic corrections to all 
orders [5] : it is therefore totally unrelated [3] to the problem of the Landau singularity discussed 
here, and we will not worry about it further. 

The minimal prescription, proposed in ref. [3], consists of defining Pll(o: s (Q 2 ),x) as an 
integral along a contour that passes to the left of the Landau pole: 



P^(a s (Q 2 ),x) = ^- T MF 
2m Jn mp - 



dNx- N lLL (a s (Q 2 ),N); < N MP < N L . (14) 



Notice that, because of the branch cut eq. 



x 



, integration to the right of the Landau pole is in 
is then free of Landau singularities. Furthermore, 



,x 



fact not possible. The function P^[ (a s (Q 

as proved in ref. [3] the series in eq. (jUJ), despite being divergent, is asymptotic to P^(a s {Q 
the difference between the minimal prescription result eq. (J14)) and the k-th order truncation of 
the divergent series eq. © is 0(a k+1 ). Interestingly, the remainder grows less than factorially 
(as (lnfc) fc for large k). 

We would like instead to tackle directly the divergent perturbative series eq. Q. The Mellin 
inversion integral can be computed explicitly: 



2m' Jn 



N+ioo 



dNx 



-N t k 



d k 1 



drj k 2m Jw-ioo 



N+ioo 



dNx~ N N~ v 



d k 



dr) k T(rj) 



In"' 



x 



where the last equality follows from the identity Jq dx x N 1 



r/=0 

+ 5(l-x) 



r?=0 



ln r '- x i 

X 



T^N^-l). Hence, 



PllMQ 2 ),x) 



9i^[-Poa s (Q 2 )] k+1 (k+1 




gij^[-(3oa s (Q 2 )] k+1 



k+l 




n=0 



k+l 



df] n T(r]) J dr] k+1 - 




In" 



-i 1 



x 



Tj=0 



j \ k+l—ri 

In In — 

x 



+6(1 -x) 



(15) 
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where in the last step we have defined A(z) = 1/T(z), and we have used the identity A( fc )(0) 
fcA^ fc_1 )(l). With straightforward manipulations we get 



PllMQ 2 ),x) 



f EhA« s (Q 2 )] n+1 

Po 



K 



ra=0 
1 



n! 



A' 



In i ^ (fc 



1. 1 1 

— [-/5 « s (Q 2 )lnln- 
n ! x 



k—n 



n + l 

K i 



5(1 



| E J A H (i: 
Po n=0 I 

+0(a 



^[-^(Q'ln-)] 
in £ x 



n+l 



+ 



-0 o a s (Q 2 )} n+] 
n + 1 



5(1 -x) 



K+l\ 



(16) 



In the limit K — > oo the terms of order af +l can be neglected, but the series is divergent. 

In the large x limit, In - = 1 — x + 0((1 
eq. (fTtjj) as 



X) 



), so, to logarithmic accuracy we may rewrite 



Pll(« s (Q 2 ),x) 



9i 

0o 



R(a s (Q 2 ),x) 



x 



(17) 



where 



K 



R(a s (Q 2 



,x 



lim ^A(")(l)[-/? a,(Q 2 (l-x))] 



n+l 



K—fOO 



n=0 



which holds up to non-logarithmic terms. Note that eq. ()18j) only follows from eq. Ijlfij) when 
X — > oo: even when is finite an infinite number of terms in eq. (|16J) is needed in order to 
reconstruct a s {Q 2 {l — x)). Equation (fTHj) can also be obtained directly by computing the inverse 
Mellin transform eq. (|15j) to logarithmic accuracy, thanks to the result of ref. [5] 



N+ioo 

— I dNx 



2ni Jn 



-N rfc+1 



k 

E 

n=0 



k + 1 

n 



A (n) (l) (ife + 1- 



n 



ln fc - n (l 



x 



(19) 



Equation (|18|) shows that the divergence of the series eqs. (J9I16)) is due to the Mellin inversion 
of IniV to all logarithmic orders: if eq. ()19|) is truncated to any finite logarithmic order, the 
resummed result in x space converges, with finite radius x < xl, as in the leading ln(l — x) case, 
eq. ()13|) . but if all logarithmic orders are included, then the series diverges, and the inclusion of 
power suppressed terms does not bring in any new divergence. 

Having understood the origin of the divergence, we can now proceed to its summation by the 
Borel method. Since we are interested in the large x limit, we neglect power-suppressed terms, 
and we use the all-log result eq. (|18|) . Namely, we take the Borel transform of the divergent 
series (JTHj) with respect to (3oot s {Q 2 {l — x)), thereby obtaining the Taylor series expansion of the 
function A(z) about z — 1: 



R(w, x) 



•E 

3=0 



A0(1) 



1 



-w 



w 



(20) 



Because A (z) is an entire function, the radius of convergence of the Borel transformed series 
eq. (j2*0j) is infinite. 
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The Borel sum of the original series is given by the inverse Borel transform 

r+oo _ w | 

R B {a s {Q 2 ),x) = - dwe /w<^a-)) — (21) 

Jo 1 (1 — w) 

However, the integrand diverges as w — > oo, because the reflection formula 

1 1 

-T(w) sm(nw) (22) 



T(l — W) 7T 

implies that A(l — w) oscillates with a factorially growing amplitude as w —>■ oo on the real axis: 
hence, the Borel sum is ill-defined. 

One may think that, alternatively, we could have performed a Borel transform with respect 
to — (3 a s (Q 2 (l — x)). In this case, we end up with 



RbMQ 2 ),x) = - / dwe*>°*W-*» — — (23) 
Jo 1 (1 + w) 

The integral now converges, because of the factorial damping provided by fn+w) as w — > oo. 
However, eq. ()23)1 diverges in the limit a s (Q 2 ) —>■ in the physical region where a s > 0. It 
is amusing to note that in the unphysical region a s < it can be easily proved that this 
Borel summation coincides with the minimal prescription. However, the physical region and 
the unphysical region of the minimal prescription manifestly cannot be analytically continued 
into each other, because in the unphysical region the cut is actually to the left of the path of 
integration in eq. (jHJ). m the physical region, instead, this modified Borel result eq. (j2SI) is 
physically unacceptable because it blows up in the perturbative limit — and in fact it is very 

w 

large even for moderate values of a s (Q 2 ), because the factor e^st^t 1 -*)) is huge before the 
damping due to the factor of T(l + w)' 1 sets in. Hence, we conclude that the result eq. (J2l!j) is 
unphysical, and we must stick with the result eq. (J2TJ). 

The presence of singularities along the path of integration in Borel inversion is a common 
occurrence in perturbative QCD, e.g. in the context of renormalons [9], and it is dealt with by 
cutting off the singularity. In our case, the singularity is as w — * oo, hence we must introduce 
an upper cutoff C to the integral. We therefore replace the divergent result eq. (}2~Tj) by 

rC | w 

RbMQ 2 ),x,C) = - dw— e /W^(i-»», (24) 

Jo 1 1 — w) 



which is convergent for all finite C . The regulated result eq. (j24|) is well defined for all x. 
Indeed, if we expand the integrand according to eq. (pjj). the series converges uniformly over the 
integration range for all finite C, so we may integrate term by term, with the result 

oo 

R B (a s (Q 2 ),x,C) = Y,^ (k) ^)[-(3ocys(Q 2 (l-x))} k+1 f k (25) 

k=0 

where 

f - 7 ( fc + 1 '/W<y a-*))) (cy^ 

fk = TikTT) (26) 
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and 7(fc, z) is the truncated Gamma function 

7(jfe+l,z) = / dwe~ w w k . (27) 
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The series eq. ({25} is convergent; however, if the cutoff C is taken to infinity, limc^oo fh = 1, 
and the original divergent series is reproduced. 

It is easy to see that the sum of the convergent series eq. (j2*5Jl is an asymptotic sum of the 
divergent series eq. (J2"T|) . Indeed, rewrite the series eq. (|23|) as 

R B (a s (Q 2 ) } x, C) = R lt (a s (Q 2 ),x) - R ht (a s (Q 2 ),x, C), (28) 

where 

oo 

Rit(a s (Q 2 ), x) = £ ^ k \l)[-p a s (Q 2 (l - x))] k+1 = R(a s (Q 2 ), x) (29) 

A:=0 

and 

(30) 

The difference between the convergent series eq. (|25|) and the first ko orders of the divergent series 
eq. (|2*T|) is equal to the sum of two terms. The first is Ru with only terms with k > ko included. 
This is of order a ko+1 . The second is the Rht, which is proportional to exp [— l/(p a s )), and 
therefore as a s — > it vanishes faster than any power of a s . The remainder of the asymptotic 
sum grows like the coefficients A^(l) of eq. (j29j). These grow less than factorially (like the 
remainder of the minimal prescription), because the Taylor expansion of the function A(z) has 
infinite convergence radius. 

The sum Rb(o:s(Q 2 ),x,C) of the series eq. ESI is regular for all < x < 1. In particular, 
it is regular at the Landau pole, where it takes a finite value which depends on C. Using the 
standard leading-order expression of a s 

<Q 2 ) = — ^ (31) 

P0 111 TT 



and the identity 

1 (k + l,z) = k\(^l-e- z J2^ zn )i (32) 
the k— th order contribution to the series eq. (|25|) is 



k x 



n=0 nl 



l -p oa , m i- x)) ri U= (--J^j e-^J^ifcin^lf^)". (33) 
In the limit x — > xl, In Q — > 0, and we get 



-/3 « s (g 2 (l - x L ))] fc+1 / fc = C fc+1 . (34) 




Figure 1: Dependence of the truncated Borel integral Rb eq. on the cutoff C for a s (Q 2 ) = 0.25 
and three values of x: below, at and above the Landau pole. 



Hence 

R B (a s (Q 2 ),x L ,C) = Y,^ k \l)jr^C k+1 = - d V A(l- V ). (35) 
t^o (k + l)\ Jo 

The same result is obtained by simply taking the limit x — > xl in eq. ()24|l. 

The dependence of the value of Rb(c(s(Q 2 ), x, C) on C is an 0(exp— — ) ambiguity in the 
definition of the asymptotic sum of the divergent series. This dependence is displayed in fig. [0 
for three values of x, below, at, and above the Landau pole, which at leading order is located at 

x L = 1 - e'Po^W). (36) 

For ease of reference, the dependence of the position xi of the Landau pole on the value of 
a s (Q 2 ) is displayed in fig. 121 Firstly, it is clear from eq. that, because 

^R B (a s (Q 2 ),x,C) = - ^ 6 -/wqV*)), (37) 

the integral has its first stationary point at C — 1, and then it is stationary at all positive integer 
values of C, where T(l — C) has simple poles. The first two stationary points are clearly visible 
in the figure. For x < xl, however, Rb(o: s (Q 2 ), x, C) has a plateau for C <> 1. The origin of this 
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0.999 — 



0.9999 — 



0.99999 — 




«.«r) 



Figure 2: Position of the leading-order Landau pole eq. as a function of a s (Q 2 ). 



i 



r(i-c) 



£ 1 

2\ 



plateau is clear from inspection of eq. (J37|) again: using eq. (|2"2"J) it is apparent that 

for C ~ 1, whereas as soon as x < xl, [3oa s {Q 2 (l — x)) < 1 in the perturbative (large Q 
c 

region, so e ^o^c^t 1 -*)) < 1. It is only when C ~ exp ^ q ^nzz^yj ^ 1 that the factorial growth 
catches up with the exponential damping. Up to this value of C the growth of _Rb(«s(Q 2 ), x, C) 
with C is negligible. This is as it should be: when x is below the Landau pole, the asymptotic 
sum of the series is essentially independent of the value of C, unless one chooses an unnaturally 
large or small value. At the Landau pole, a dependence on the value of C appears, due to the 

c 

fact that if x <* xl the exponential 0(exp(l/a s )) prefactor e Aj^cs 2 ^-*)) ^ 1, and the concept 
of asymptotic sum starts losing its meaning. A minimal choice of C is C — 1, so Rb is stationary 
for all x, and in the beginning of the plateau for x < xl- 
The physical meaning of C becomes clear by rewriting 



' l3 a s (Q 2 (l-x)) 



A 2 



x 



c 



(38) 



In other words, Rht(a s (Q 2 ), x, C) eq. (fHU|) is a higher twist contribution: the original divergent 
perturbative series eq. (|29jl has been made convergent by the inclusion of the higher twist series 
eq. (J30J). Of course, eq. (|2*5jl implies that this twist series must necessarily also be divergent. 
Integer values of C correspond to even twists, and in particular if C — 1, Rht{c( s (Q 2 ),x,C) is a 
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standard twist-4 contribution, namely, the first subleading twist. The choice C = 1 is minimal 
in that it corresponds to regulating the Borel summation through the first subleading twist. 

Equation (|3*5Jl implies that these higher twist terms are suppressed by powers of 4a, but 
enhanced by powers of j^-. The Landau pole is the point where the parameter of the twist 
expansion is equal to one, i.e., leading twist and higher twist terms are of comparable size. 
However, despite this enhancement, as long as we choose C < 1 the Borel sum eq. ()24|) remains 
integrable at x = 1. Indeed, with C = 1 we have 

rl | w 

{l-x)R B {a s {Q 2 ),x,l) = - dw— -e ^W^> (1 - x) l ~ w , (39) 

Jo T(l — w) 

which vanishes as x — > 1 thanks to the fact that < w < 1. It follows that Pll eq. (jl7|) 
with R given by R-q eq. (|24|) acts as a conventional + distribution, and in particular the integral 
Jo 1 PhhictsiQ 2 ), x)f(x)dx is finite if the test function f(x) is regular at x — 1. 

When Rm is viewed as a genuine higher twist term, the prefactor ^ in eq. (jHOl) comes 
from the Wilson expansion, and not from a factor of exp — , so the higher twist term is just of 
0(a s ). However, this term matches an ambiguity in the leading twist, which does not appear 
at any finite order in the expansion of the leading twist term itself but only in its asymptotic 
Borel resummation. Equivalently, the higher twist contribution removes the cutoff ambiguity 
introduced by the need to treat this divergence. The situation is thus akin to the customary 
case of renormalons, where similarly the ambiguity introduced by the need to make the Borel 
inversion well-defined is cured by the inclusion of higher twist terms. Henceforth, we will take 
our result to be given by the Borel sum eq. ()24|) with C = 1. 

We now wish to compare this result to the minimal prescription. We do so by defining 

RmpMQ 2 ),x) = (l-x)P MP (a s (Q 2 ),x), (40) 

where Pmp(&s{Q 2 ), x ) is given by eq. (fT4^ . It is clear that the Borel asymptotic sum eq. (j^j) and 
the minimal prescription asymptotic sum eq. (J4U|) of the divergent series eq. (|21|) cannot coincide 
in general, because the former depends on C and the latter doesn't. Their x-dependence when 
C — 1 is compared in fig. El where we also show the value of R obtained by truncating the 
perturbative expansion of the divergent series eq. (J21j) . This is defined by including in R the 
contribution from the first K terms in the expansion in powers of a s (Q 2 ) eq. (jHJ), where K is 
defined as the value where the k th term in the sum in eq. Q, ctfc, starts growing, i.e. as the value 
K such that 

\a k + a k -i\ > \a k - 2 + a fc _ 3 | for k > K . (41) 

The definition of the optimal truncation point has some degree of arbitrariness; we have checked 
that with our choice, eq. (j4*T|) . the truncated sum is closer to the asymptotic sum than it would 
be by simply requiring that \a k \ > \a k -i\ for k > K . Note that we have defined the truncation 
in terms of the expansion eq. in powers of a s (Q 2 ), and not in terms of that eq. (|21j) in powers 
of a s {Q 2 {\ — x)), in order for the truncated result to be well-defined also at the Landau pole xl 
where a s {Q 2 {l — x)) blows up. For comparison, we also include in fig. El the leading ln(l — x) 
result eq. ()13|). and the large x form of the unresummed result, which is simply given by 

RloMQ 2 ),x) = -(3 a s (Q 2 ) + O[(l-x)]. (42) 
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Figure 3: Various determination of the x-space resummed result for a s (Q 2 ) = 0.25. LLx denotes 
the leading ln(l — x) result of eq. f|13|) . while B, MP and pert denote three different determinations of 
the divergent leading IniV series eq. (|18|). respectively through Borel summation eq. (|24() with C = 1, 
the minimal prescription eq. (|14fl and the asymptotic truncation of the perturbative expansion at K$ 
eq. (|41l) . The large-x (constant) leading order result eq. ()42|) is also shown for comparison. 



Below the Landau pole the Borel and minimal resummation prescriptions are close to the 
truncated perturbative result and thus close to each other, as one expects of an asymptotic sum. 
Note that the leading ln(l — x) result eq. (|T3*|) is reasonably close to these results but never quite 
on top of them: at not so large x, where it reduces to the leading-order result eq. (jUJ), it differs 
from them because of subleading terms — this is the region were the large x resummation is not 
very useful. As one enters the large x region, however, the leading ln(l — x) result eq. (fT^j) is 
contaminated by the Landau pole where it blows up. So the leading ln(l — x) result turns out 
to be of limited usefulness, because there is no region of x where it is applicable. At and above 
the Landau pole the Borel prescription and the MP prescription start deviating: in this region 
the higher twist contributions which stabilize the Borel sum are of the same order as the leading 
twist. On the other hand, at and above the Landau pole the series diverges very fast and its 
asymptotic sum looses meaningfulness. Hence, comparison of the two prescriptions (Borel and 
minimal) gives us an estimate of the size of nonperturbative effects: when the two prescriptions 
start departing from each other, nonperturbative effects become important. Indeed, these two 
prescriptions bracket the truncated perturbative expansion, which oscillates between them as 



II 



the order of the truncation (the value of K ) varies as a function of x. 

In summary, we have traced the origin of the divergence of the perturbative expansion of 
soft gluon resummation, and we have shown that it may be treated by Borel resummation 
stabilized by higher twist terms. The result that we found is close to the widely adopted minimal 
prescription, but it deviates from it when nonperturbative corrections become important, namely 
at the Landau pole. All our computations were presented in the case of threshold resummation 
(such as e.g. DIS at large x) at the leading logarithmic level. The extension to all logarithmic 
orders and to pr resummation will be discussed elsewhere. Our result is useful for practical 
calculations in that it does not require the numerical evaluation of a Mellin inversion integral. 
Furthermore, the availability of more resummation methods that differ in the nonperturbative 
region is useful in order to assess the reliability of perturbative resummed results. 
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